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Abstract. The formation of the UV OH spectral lines has been investigated for a range of stellar parameters in the light of 
3D hydrodynamical model atmospheres. The low atmospheric temperatures encountered at low metallicities compared with the 
radiative equilibrium values enforced in classical ID hydrostatic model atmospheres have a profound impact on the OH line 
strengths. As a consequence, the derived O abundances using 3D models are found to be systematically lower by more than 
0.6 dex at [Fe/H] = —3.0 compared with previous ID analyses, casting doubts on the recent claims for a monotonic increase 
in [O/Fe] towards lower metallicities. In fact, taken at face value the resulting 3D LTE trend is in rough agreement with the 
. conventional [O/Fe] plateau. Caution must, however, be exercised in view of the remaining assumptions in the 3D calculations. 

We have verified that the stellar parameters remain essentially unchanged with 3D model atmospheres provided that the infrared 
flux method (AT c ff 20 K), Hipparcos parallaxes (Alogp 0.05) and Fell lines (A[Fe/H] 0.1 dex) are utilised, leaving 
the 3D O abundances from OH lines largely intact (A[0/H] £ 0.05 dex). Greater concern stems from possible departures 
from LTE in both the line formation and the molecular equilibrium, which, if present, would increase the derived O abundances 
again. Non-LTE line formation calculations with ID model atmospheres suggest no significant steepening of the [O/Fe] trend 
even if the abundance corrections amount to about 0.2 dex for all investigated stellar parameters. We note, however, that the 
3D case may not necessarily be as metallicity-independent. The apparent lack of laboratory or theoretical rate coefficients at 
the relevant temperatures for the involved molecular reactions unfortunately prevents a quantitative discussion on the possible 
effects of non-equilibrium chemistry. 
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& ' 1 . Introduction elements Li, Be and B are produced through cosmic ray spal- 

lation of C, N and, most importantly, O with protons and a- 
Oxygen is the third most abundant element in the Universe after particles xhuSi a proper unders tanding of the oxygen abun- 
hydrogen and helium. Besides being relatively common, oxy- dancefj of stars of different metallicities is required in order to 
gen has also attracted a great deal of attention due to its special interpret ±s evolution of the light elements, 
role in galactic and stellar evolution. Since oxygen is predom- 

. . • t u ii c - Following the first indication of an oxygen over-abundance 

inantly produced in connection with core collapses of massive ° JD 

,. TT OKT TT s ,., • , £ , , relative to iron ([O/Fe] > 0) in metal-poor stars (Conti et al. 

stars (type II supernovae, SNe II) while iron is also forged dur- * 

. , , „ „ ... , . ,, T 1967), a great number of studies have been devoted to quan- 

mg accretion-induced collapses of white dwarfs (type la super- , . 

ox , j \ i i i . tify this enhancement. Although all agree on its existence, the 

novae, SNe la), the abundance ratio of these two elements at J e & 

j.™ u • ,l c amount of the over-abundance is hotly contested. Various oxy- 

different cosmic epochs give important insight to the forma- . J 

, , .■ c ^ , „ . . . gen diagnostics in different types of stars have been applied 

tion and evolution of the Galaxy, as well as constraining the & & ->r rr 

, ■ c rr,, . t ■ , but with disparate results. The forbidden [O I] lines at 630.0 

physics of supernovae. The amount of oxygen in stars can also r . . L J 



and 636.3 nm in metal-poor giants suggest a nearly flat plateau 
at [O/Fe] - 0.4 for [Fe/H] <, -1.0 (e.g. Barbuy 1988; Sneden 
et al. 1991) while the O I niplet at 777 nm in metal-poor dwarfs 



significantly influence the nuclear energy production and the 

opacities, which affects stellar evolution. As a consequence, the 

dating of globular clusters depend on the stellar oxygen con- , 

, , r^v/r-nlTlr and subgiants tend to imply systematically higher values, of - 

tent: an increase of the oxygen over-abundance [O/Fe] from ° r j j j a 

, „ „ . , 1 „ . , u i r.uu ten with a monotonic increase towards lower metallicities (e.g. 

+0.3 to +1.0 introduces an about 2 Gyr lower age of the old- ° 

. i . /\r j n on 11 onnn a jj-.' n i- u. Abia & Rebolo 1989; Israelian et al. 1998, 2001 ; Boesgaard et 

est clusters (VandenBerg & Bell 2001). Additionally, the light , ,„„„ „ ' „„„„, b 

5 3 ' 5 al. 1999; Carretta et al. 2000). 



Send offprint requests to: e-mail: martin@astro.uu.se All oxygen criteria have their pros and cons, which influ- 

1 The abundance ratios are defined by the customary ence the conclusions. The forbidden lines are immune to de- 
[X/Fe]=log(A'x/A r Fe)* — log(7Yx/-Wpe)Q partures from local thermodynamic equilibrium (LTE) (cf. dis- 
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cussion in Kiselman 2001) but the lines are very weak at low 
metallicities, in particular in unevolved stars where the feature 
becomes undetectable for [Fe/H] ^ —2.0 (Nissen et al. 2001). 
Furthermore, concerns regarding the primordial nature of the 
oxygen in field giants have been voiced as mixing of nuclear- 
processed material may pollute the surface, as evident in many 
globular cluster giants (e.g. Langer et al. 1997). The O I triplet 
on the other hand is more easily discerned in metal-poor dwarfs 
but also more susceptible to departures from LTE (cf. Kiselman 
2001 and references therein) and inhomogeneities introduced 
by stellar granulation (Kiselman & Nordlund 1995) than the 
[O I] lines. The high excitation potential of the O I lines also 
make them vulnerable to errors in the effective temperature 
(Toff)- Additionally, an often overlooked source of confusion 
is the stellar Fe abundances, which should be self-consistently 
computed in order to obtain reliable [O/Fe] ratios. Finally, even 
the question of the absolute solar oxygen abundance is still 
not settled, leaving the fundamental reference point for stellar 
[O/Fe] ratios insecurely anchored. 

An attractive alternative to the O I and [O I] lines are pro- 
vided by the molecular OH electronic transitions in the UV. 
(Bessell et al. 1984, 1991; Nissen et al. 1994). In fact, at 
[Fe/H] ~ — 4 no substitute for the OH lines as the prime oxygen 
diagnostic is available. Therefore the OH spectral line forma- 
tion must be properly understood when attempting to probe the 
earliest epochs of the Galactic evolution for example to iden- 
tify the nucleosynthetic fingerprints of the elusive Population 
III stars (e.g. Karlsson & Gustafsson 2001). Recently, Israelian 
et al. (1998, 2001) and Boesgaard et al. (1999) have analysed 
the UV OH (A 2 £ - X 2 n) lines in stars down to [Fe/H] ~ -3.3 
with high S/N and resolution spectra and have found a lin- 
ear trend in [O/Fe] vs [Fe/H] with a slope of about —0.4^, 
in stark conflict with the [O I] results. If confirmed, these re- 
sults would have far-reaching consequences, as outlined above. 
The situation is complicated, however, by recent studies of the 
OH vibrational-rotational lines in the infrared (IR), which sug- 
gest a nearly flat [O/Fe] in agreement with the [O I] findings 
(Balachandran et al. 2001; Melendez et al. 2001). 

As all previous investigations have been based on ID hy- 
drostatic model atmospheres, one may worry about possible 
systematic errors introduced by the inherent assumptions of 
the analyses. Recently the first 3D hydrodynamical model at- 
mospheres of metal-poor stars have been constructed (Asplund 
et al. 1999a), which have very different temperature structures 



When restricting to the OH-based results and taking the errors in 
both [Fe/H] and [O/Fe] into account, the Boesgaard et al. (1999) data 
gives a slope of -0.40 ± 0.04 and -0.36 ± 0.04 for the King (1993) 
and Carney (1983) T c ff-scales, respectively. Similarly, the published 
results of Israelian et al. (1998) are consistent with a slope of — 0.38± 
0.07, which becomes —0.37 ± 0.06 when also including the new OH 
data in Israelian et al. (2001). For simplicity we here adopt a slope of 
—0.40 but emphasize that the exact value only has a marginal effect 
on our results. It should be noted that the values given in the original 
references are slightly smaller due to the combination of OH and O I 
results (Boesgaard et al.) and restriction to stars with [Fe/H] < —1.0 
(Israelian et al.). When also considering the non-LTE results for Fe I 
by Thevenin & Idiart (1999), the slopes decrease by about 0.06 (King 
2000; Israelian et al. 2001), but cf. discussion in Sect. 5.4. 



compared with classical ID model atmospheres. As a con- 
sequence, analyses of temperature sensitive spectral features 
can be suspected to be systematically in error if relying on 
ID model atmospheres. In particular, Asplund et al. (1999a) 
cautioned that oxygen abundances of metal-poor stars derived 
from ID studies of OH lines may be strongly overestimated 
due to the extreme temperature sensitivity of molecule forma- 
tion. The aim of the present paper is to quantify this suspicion 
in terms of O abundances and investigate other possible sys- 
tematic errors which may hamper the 3D analysis of OH lines. 
Preliminary results have been presented in Asplund (2001). 



2. 3D hydrodynamical model atmospheres 

Realistic ab-initio 3D, time-dependent simulations of stel- 
lar surface convection form the foundation for the present 
study. The same incompressible radiative hydrodynamical code 
which previously has been applied successfully to studies of so- 
lar (e.g. Stein & Nordlund 1998; Asplund et al. 2000a,b) and 
stellar granulation (e.g. Asplund et al. 1999a, Allende Prieto 
et al. 2001; Asplund et al., in preparation) has here been used 
to construct sequences of 3D model atmospheres with vary- 
ing stellar parameters. The equations of mass, momentum and 
energy conservation together with the simultaneous treatment 
of the 3D radiative transfer equation have been solved on a 
Eulerian mesh with 100 x 100 x 82 gridpoints. The physical di- 
mensions of the numerical grid were sufficiently large to cover 
many 10) granules simultaneously. The depth scales have 
been optimized to provide the best resolution where it is most 
needed, i.e. in those layers with the steepest gradients in terms 
of dT/dz and d 2 T/dz 2 , which for the solar-type stars occurs 
around the visible surface. 

Special care has been exercised to include the most appro- 
priate input physics. In particular, state-of-the-art equation-of- 
state (Mihalas et al. 1988), which includes the effects of ioniza- 
tion, excitation and dissociation of the most important atoms 
and molecules, and relevant continuous (Gustafsson et al. 1975 
with subsequent updates) and line (Kurucz 1993) opacities 
have been employed. During the convection simulations, the 
3D radiative transfer is solved for in total eight inclined rays un- 
der the simplifying assumptions of LTE (S v = B v ) and group- 
ing of the opacities into four bins (Nordlund 1982). At regu- 
lar intervals during the simulations, the accuracy of the opacity 
binning technique is verified by solving the full monochromatic 
radiative transfer (about 2700 wavelength points) in the 1 .5D 
approximation, i.e. treating each vertical column as a sepa- 
rate ID model atmosphere and ignoring all horizontal radiative 
transfer effects. Further details on the numerical procedures of 
the simulations may be found in Stein & Nordlund (1998). 

For the present purpose, two sequences of 3D model at- 
mospheres have been constructed. The first series of mod- 
els (here: the solar sequence) correspond to the Sun (T c g ~ 
5800K [J log g — 4.44 [cgs]) but with a range of metallicities 



Since the entropy of the inflowing gas at the lower boundary has 
replaced T c g as an independent input parameter in 3D convection sim- 
ulations, the resulting T c s varies slightly in time due to the evolution 
of individual granules. In order to obtain a specific temporally aver- 
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Table 1. Details of the 3D hydrodynamical model atmospheres 



< Toff > a 


log g 


[Fe/H] 


x,y,z-dimensions 


time b 


[K] 


[cgs] 




[Mm] 


[min] 


5767 ±21 


4.44 


±0.0 


6.0x6.0x3.7 


50 


5822 ± 11 


4.44 


-1.0 


6.0x6.0x3.7 


50 


5837 ± 8 


4.44 


-2.0 


6.0x6.0x3.7 


40 


5890 ± 8 


4.44 


-3.0 


6.0x6.0x3.7 


40 


6191 ±37 


4.04 


±0.0 


21.4x21.4x8.7 


60 


6180 ± 20 


4.04 


-1.0 


21.4x21.4x8.7 


60 


6178 ± 17 


4.04 


-2.0 


21.4x21.4x8.7 


60 


6205 ± 15 


4.04 


-3.0 


21.4x21.4x8.7 


60 



a The temporal average and standard deviation of the emergent T c a 
of the different convection simulations 

b The time coverage of the part of the simulation used for the spec- 
tral line calculations, while the full convection simulations are 
much longer extending over several convective turn-over time- 
scales 



([Fe/H] = 0.0, -1.0, -2.0 and -3.0), while the second suite 
(here: the turn-off sequence) correspond to typical turn-off stars 
(T c ff ~ 6200 K, log g = 4.04 [cgs]), again with varying metal 
content ([Fe/H] = 0.0, -1.0, -2.0 and -3.0). Some details 
of the simulations are given in Table [j]. The individual ele- 
mental abundances have been taken from Grevesse & Sauval 
(1998) scaled appropriately to the relevant [Fe/H] with no a- 
element enhancements. As long as the comparison is strictly 
differential between the 3D and ID model atmospheres as in 
our case the omission of a-element enhancements will be neg- 
ligible for the spectral line formation. The initial snapshots for 
the simulations were taken from simulations of lower numer- 
ical resolution (50 x 50 x 82), which had been run sufficiently 
long times to allow thermal relaxation and a statistically steady 
state at the wanted T c ff to be established. The initial snapshots 
for these lower resolution runs in turn were obtained from a 
previous solar simulation (Stein & Nordlund 1998) scaled ap- 
propriately to the new stellar parameters using the experience 
from ID hydrostatic stellar models and the entropy variations 
in 2D hydrodynamical model atmospheres (Ludwig et al. 1999; 
Freytag et al. 1999). Classical ID, hydrostatic MARCS model 
atmospheres (Asplund et al. 1997) with identical input param- 
eters and chemical compositions as the 3D simulations have 
been constructed to allow a differential comparison in terms 
of spectral line formation. Whether 3D models indeed should 
be assigned the same T e s as a cor responding ID model atmo- 
sphere is investigated in Sect. 5.2 when comparing the predic- 
tions for the infrared flux method (IRFM) with the two types of 
models. 



aged value for the emergent T c s , a careful and very time-consuming 
fine-tuning of the inflowing entropy would be required. Since we are 
only interested in a differential comparison between 3D and ID, we 
have not attempted to obtain exactly the solar T c ff for the individual 
simulations and instead settled for values in reasonable proximity of 
the targeted T c s , cf. Table |l| Naturally, the comparison ID model at- 
mospheres have the same T c g as the final 3D average to isolate the 
granulation effects from differences in input parameters. 



Although the details of the convection properties of the 
present and other simulations for late-type stars will be de- 
scribed elsewhere (Asplund et al., in preparation), the most im- 
portant metallicity effects on the resulting photospheric struc- 
tures are briefly discussed here in order to understand the im- 
pact of the 3D models on the spectral line formation of the OH 
lines. While the temperature remains close to the radiative equi- 
librium value at solar metallicities and mild metal-deficiencies 
([Fe/H]> —1.0), the temperature in the outer layers depart sig- 
nificantly from it at lower metallicities (Asplund et al. 1999a). 
The temperature in this optically thin region is determined from 
a competition between adiabatic cooling and radiative heating. 
The latter mechanism arises when continuum photons released 
at deeper layers are reabsorbed in spectral lines. At solar metal- 
licities the abundant spectral lines succeed in providing suf- 
ficient radiative heating to balance the adiabatic cooling and 
keeping the average temperature close to the radiative equi- 
librium expectation. At progressively lower metallicities, the 
available spectral lines become fewer and weaker, which al- 
lows cooling to dominate more. As a consequence, balance is 
not restored unless the temperature is much below the radia- 
tive equilibrium value. Thus, in hydrodynamical model atmo- 
spheres spectral lines have the opposite effect compared to in 
hydrostatic model atmospheres which enforce radiative equi- 
librium: spectral lines cause surface heating. The effects of the 
lower temperatures are first visible at the outermost layers but 
move towards deeper layers at progressively lower metallici- 
ties. 

When viewed as a function of optical depth instead of ge- 
ometrical height, the low surface temperatures are less pro- 
nounced since the temperature sensitivity of the continuous 
opacities partly hides the effect (Asplund et al. 1999a); lower 
temperatures also imply smaller opacities and optical depths. 
Nevertheless the low temperatures reach typical line-forming 
regions at low metallicities. Furthermore, the cool outer lay- 
ers shift the whole line formation outwards for temperature 
sensitive features such as molecular lines. For example, at 
log T500 = — 3 the average temperature difference between 3D 
and ID models can reach 1000K and thus have a profound ef- 
fect on the lines sensitive to those layers. 

Differences in the gas and electron pressures could also in- 
fluence spectral line formation. With the exception of some ad- 
ditional contribution from turbulent pressure in the convective 
overshoot region around the visual surface, the total pressure is 
well approximated by the gas pressure in the photosphere for 
late-type dwarfs. Since even in the 3D hydrodynamical models 
the photosphere is typically not too far from hydrostatic equi- 
librium, the pressure scale height Hp = — (dlnP/dr) -1 
P gas /gp oc T and thus the resulting 3D pressure structure tend 
to be lower than the corresponding ID value at a given geo- 
metric height at low metallicities. Naturally, the lower temper- 
atures have an even greater impact on the electron pressure. 
The lower gas and electron pressures in metal-poor 3D model 
atmospheres will affect mainly lines that are considered gravity 
sensitive. 
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3. Spectral line formation in 3D model 
atmospheres 

The OH line calculations with the 3D hydrodynamical model 
atmospheres follow the same procedure as in other recent in- 
vestigations of the influence of granulation on stellar spec- 
troscopy (Asplund et al. 1999a,b, 2000a,b,c; Asplund 2000a,b, 
2001; Asplund & Carlsson 2001; Allende Prieto et al. 2001; 
Nissen et al. 2000, 2001; Primas et al. 2000). The 3D con- 
vection simulations, which extend down to the essentially adi- 
abatic layers well below the visible photosphere, were inter- 
polated to a finer vertical depth-scale although with the same 
number of depth-points to improve the numerical accuracy in 
the spectral synthesis. The average continuum optical depth 
range typically beyond logr 50 o ^ 2 and above logr 50 o ^ 
—5, in order to minimize the influence from the artificial top 
and bottom boundaries. Prior to the line transfer calculations 
the horizontal resolution was decreased from 100x100x82 
to 50 x 50 x 82 to ease the computational burden; various tests 
ensured that the procedure did not introduce any differences 
in the spatially and temporally averaged line profiles. From 
the full convection simulations, which cover several convective 
turn-over time-scales, representative sequences of typically one 
hour stellar time with snapshots every 30 s were selected for 
the spectral syntheses. In terms of derived oxygen abundances, 
the number of snapshots were sufficient to provide statistically 
significant results (Alogeo < 0.01 dex), as verified by test cal- 
culations with shorter time sequences. The mean T c ff 's for the 
shorter simulation sequences are given in Table [l]. 

The spectral line formation in 3D model atmospheres was 
performed under the assumption of LTE. Thus, the OH number 
densities were computed from instantaneous molecular equilib- 
rium and Saha ionization and Boltzmann excitation balances. 
The line source function was approximated with the Planck 
function (S u = B u ). The applicability of thes e st rong as- 
sumptions is further investigated in Sects. 5^5 and 5J5 The OH 
molecular equilibrium was computed with scaled solar abun- 
dances with the exception of the oxygen abundance. The re- 
moval of available oxygen atoms due to CO formation was 
taken into account but the effect was found to be negligible. 
Flux profiles were computed for two typical OH UV lines 
(313.9 and 316.7 nm) from solving the radiative transfer for 
in total 17 inclined rays (four /i-angles and four </?-angles plus 
the vertical /i = 1.0). The line transition data for the OH lines 
was taken from Israelian et al. (1998); we emphasize that in a 
differential 3D- ID comparison such as ours the exact choices 
of for example the ^/-values are not important. In addition, 
we have included six of the Fe II lines used by Nissen et al. 
(2001) in their study of [O/Fe] from the [Ol] line to quantify 
the corresponding impact on Fe abundances. The background 
continuous opacities were calculated using the Uppsala opac- 
ity package (Gustafsson et al. 1975 with subsequent updates). 
Since the Doppler shifts introduced by the convective velocity 
field are fully accounted for, no microturbulence or macrotur- 
bulence parameters enter the 3D line calculations (Asplund et 
al. 2000b,c). It is noteworthy that none of the various free pa- 
rameters hampering ID analyses (e.g. mixing length parame- 
ters, micro- and macroturbulence) are necessary with 3D hy- 



drodynamical model atmospheres. For the ID calculations, a 
microturbulence of l.Okms -1 has been assumed for the solar 
and turn-off sequences; the choice of microturbulence, how- 
ever, is only of some significance for the solar metallicity mod- 
els while unimportant at low metallicities. 

The main advantage of limiting the comparison to a strictly 
differential study of the 3D and ID predictions is that the un- 
certainties in e.g. the absolute transition probabilities, equiva- 
lent widths, blends, continuum placement, missing UV opac- 
ities, r c ff-calibration and the solar oxygen abundance can be 
avoided. As a consequence we are not able to determine ac- 
curate absolute stellar oxygen abundances. Instead we only at- 
tempt to address the question of whether there are systematic 
errors in analyses of the UV OH lines when relying on classical 
ID model atmospheres. A star-by-star analysis of high quality 
OH observations is left for a future paper. 



4. 3D LTE oxygen abundances derived from UV 
OH lines 

Molecule formation is extremely temperature sensitive (as long 
as the molecule in question is a trace element for its atomic con- 
stituents, which is the case for OH in the stars studied here). For 
conditions typical of the line-forming regions in the Sun, the 
LTE number density of OH, Non, is proportional to w 7 1 - 12 
as estimated from our adopted equation-of-state and molecu- 
lar balance routines. This high degree of non-linearity makes 
molecular lines very susceptible to errors in the adopted tem- 
perature structures of the model atmospheres. In particular, the 
very different temperature structures of typical convective up- 
and downflows compared with classical ID model atmospheres 
(i.e. up- and downflows can not be represented by two theoret- 
ical hydrostatic model atmospheres with different T ff, since 
granules have a much steeper temperature gradient than inter- 
granular lanes), can be suspected to significantly influence the 
molecular number densities and line formation. 

Under the assumption of LTE adopted here, Aoh depend 
only on the instantaneous local temperature. In general the 
temperature contrast reverses in the convective overshoot re- 
gion some distance above the continuum forming layers, i.e. 
gas above the warm, upflowing granules tend to be cooler than 
average (e.g. Stein & Nordlund 1998). In particular in metal- 
poor stars the horizontal temperature contrast is very large due 
to the weak coupling between the gas and the radiation field 
(Asplund et al., in preparation). In general the temperatures at 
the high atmospheric layers are very low but occasionally the 
compression from converging gas flows or shocks can rise the 
temperature to the radiative equilibrium value or even above 
it. As a consequence, the UV OH LTE line strengths across 
the stellar granulation are expected to be stronger in the gran- 
ules but with pronounced scatter when viewed nearly face-on, 
as confirmed by the 3D line calculations (Fig. [j]). The spatially 
and temporally averaged OH flux profiles are therefore strongly 
biased towards upflowing regions due to their higher contin- 
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Fig. 1. The contribution to the equivalent width, W\(x, y) • 
^ ont (x,y)/(/ cont ), of the OH 31 3.9 nm line across the stellar 
surface forthe [Fe/H]= +0.0 (Upperpanel) and [Fe/H]= -3.0 
(Lower panel) models in the solar sequence as a function of the 
relative continuum intensity /™ nt (x, y)/(/ cont ) at 313.9nm. 
The spatially resolved intensity profiles have been computed 
for /i = 1.0. The horizontal and vertical dashed lines denote 
the mean equivalent widths and continuum intensities with 3D 
model atmospheres, while the solid lines mark the correspond- 
ing ID predictions using the same O abundances as in 3D. Note 
that the here shown 3D calculations only refer to one snap- 
shot each for the two models and not the whole simulation se- 
quences 



uum intensities, steeper temperature gradients and larger area 
coverage [| 

Typical resulting 3D and ID line profiles are shown in Fig. 
||, which clearly illustrates the influence of 3D models on the 
derived O abundances from UV OH lines at low metallicities. 
The much lower temperatures and larger concentration of OH 
molecules in the 3D model atmospheres of metal-poor stars re- 
sult in far stronger OH lines than with ID model atmospheres 
while the difference is much less accentuated at solar metallici- 



Note that the opposite is true for the OH pure rotational lines in 
the IR since these lines and their surrounding continuum are formed at 
significantly higher layers where the temperature contrast is reversed 
relative to the visual continuum (cf. Figs. 1 and 7 in Kiselman & 
Nordlund 1995) 




313.915 313.920 
Wavelength [nm] 



Fig. 2. The predicted profile for the OH313.9nm line using 
the 3D hydrodynamical model atmosphere for Teff=5800K, 
logg=4.44 [cgs], [Fe/H]=-3.0 with an adopted O abundance of 
log eo = 6.4 (solid line). In comparison the corresponding ID 
profile with the same abundance is much weaker (upper dashed 
line). Only with logeo = 7.05 (lower dashed line) is the ID 
profile comparable in strength to the 3D profile. In no cases 
have extra broadening in the form of rotation or macroturbu- 
lence been applied 

Table 2. Comparison of the oxygen LTE abundances derived 
with ID hydrostatic and 3D hydrodynamical model atmo- 
spheres. The 3D O abundances are those which reproduce the 
equivalent widths computed using a ID model atmosphere, a 
microturbulence of £turb = 1.0 km s -1 and the oxygen abun- 
dances given in the fourth column 



Toff 




[Fe/H] 


logeo, id 


log eo,3D 


log eo,3D 


[K] 


[cgs] 






313.9 nm a 


316.7 nm b 


5767 


4.44 


+0.0 


8.90 


8.88 


8.89 


5822 


4.44 


-1.0 


8.30 


8.03 


8.07 


5837 


4.44 


-2.0 


7.70 


7.19 


7.25 


5890 


4.44 


-3.0 


7.10 


6.44 


6.52 


6191 


4.04 


+0.0 


8.90 


8.80 


8.82 


6180 


4.04 


-1.0 


8.30 


7.97 


8.00 


6178 


4.04 


-2.0 


7.70 


6.90 


7.01 


6205 


4.04 


-3.0 


7.10 


6.08 


6.22 


HM, int. c 


4.44 


+0.0 


8.90 


8.69 


8.72 


HM, flux c 


4.44 


+0.0 


8.90 


8.64 


8.69 



log gf = -1.76, x CX c = 0.76 eV (Israelian et al. 1998) 
logff/ = -1.69, Xoxc = 1.11 eV (Israelian et al. 1998) 
Using the Holweger-Miiller (1974) ID semi-empirical model at- 
mosphere for the solar case instead of a MARCS solar model, ei- 
ther using intensity or flux line profiles 



ties, in accordance with the behaviour of the temperature struc- 
tures. Table || lists the O abundances of the two OH lines for 
the two types of model atmospheres. The 3D O abundances are 
those which reproduce the equivalent widths computed using 
ID model atmospheres and an adopted [O/Fe] = — 0.40-[Fe/H] 
trend (Israelian et al. 1998, 2001; Boesgaard et al. 1999). A 
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less steep [O/Fe] trend would have resulted in slightly smaller 
granulation effects due to the shifting of the line-formation re- 
gion inwards where the low temperatures in the 3D models 
are less pronounced. Similarly, the impact of the 3D models 
is greater for the OH 313. 9nm line than for the OH 316.7 nm 
line in accordance with the lower excitation potential and larger 
line strength of the former. The in general larger granulation 
corrections for the hotter models are in agreement with their 
larger temperature differences and the high non-linearity of the 
OH line formation. 

For completeness we have carried out an identical 3D- ID 
comparison for six Fe II lines. The resulting granulation correc- 
tions are listed in Table [}| As expected (Asplund et al. 1999a), 
the impact of 3D model atmospheres are relatively minor on 
the Fe II lines since they are formed in deep atmospheric layers 
and therefore not sensitive to the low temperatures encountered 
in the upper layers in the 3D models. In terms of Fe abundance, 
the difference between the ID and 3D predictions amounts to 
< 0.1 dex. 

The strong metallicity dependence of the abundance cor- 
rections has a profound impact on the use of UV OH lines as 
O diagnostic in metal-poor stars and cast serious doubts on re- 
cent claimed linear trends in [O/Fe] towards lower metallic - 
ities based on ID LTE analyses (Israelian et al. 1998, 2001; 
Boesgaard et al. 1999). In Fig. ||]the average granulation correc- 
tions of the two investigated OH lines for the solar and turn-off 
sequences are shown, while Fig. [| presents the final mean cor- 
rections of the different T c g simulations [J. The correction of the 
ID LTE result due to the use of 3D model atmospheres is here 
defined as relative to the solar calibration, i.e. the abundance 
differences in Table |^ are subtracted with the corresponding ef- 
fect for the Sun since [O/Fe] ratios are studied. This procedure 
is similar to the common practice of determining astrophysical 
^/-values from the Sun, as was also the basis for the works by 
Nissen et al. (1994), Israelian et al. (1998, 2001) and Boesgaard 
et al. (1999). For the Sun, two possible calibrations are pos- 
sible using either the Holweger-Miiller (1974) semi-empirical 
model atmosphere or the theoretical MARCS model as the ID 
representation of the solar photosphere, yielding two possible 
3D trends in Figs. || and [|. The difference between the two 
trends therefore simply reflects the difference in derived solar O 
abundances when relying on the two types of ID model atmo- 
spheres. It should be noted that the granulation corrections for 
Fe II lines presented in Tableware not included in Figs. ^JandQ, 
since the adopted [O/Fe] trend with metallicity originates from 
analyses of Fel lines (Boesgaard et al. 1999; Israelian et al. 
1998, 2001, cf. discussion in Sect. 0). According to Table [| 
the inclusion of the Fe II results would bring down the 3D LTE 
[O/Fe] results further by ^ 0.1 dex at the lowest metallicities. 

Taken at face value, with the claimed ID [O/Fe] results and 
the here presented granulation corrections, the emerging trend 
with metallicity is in fact roughly consistent with the since long 
advocated [O/Fe] plateau from [Ol] lines, as seen in Fig. 
However, we caution that such a conclusion is likely prema- 



Table 3. Comparison of the iron LTE abundances derived with 
ID hydrostatic and 3D hydrodynamical model atmospheres 
from Fe II lines. The 3D Fe abundances are those which re- 
produce the equivalent widths computed using a ID model at- 
mosphere, a microturbulence of £tmb = l.Okms -1 and the Fe 
abundances given in the fourth column 



[K] 


[cgs] 


[Fe/H] 


log £p c ,lD 


log <Tc,3D 


Alogep c 


5767 


4.44 


+0.0 


7.50 


7.48 


-0.02 


5822 


4.44 


-1.0 


6.50 


6.54 


+0.04 


5837 


4.44 


-2.0 


5.50 


5.60 


+0.10 


5890 


4.44 


-3.0 


4.50 


4.59 


+0.09 


6191 


4.04 


+0.0 


7.50 


7.49 


-0.01 


6180 


4.04 


-1.0 


6.50 


6.52 


+0.02 


6178 


4.04 


-2.0 


5.50 


5.56 


+0.06 


6205 


4.04 


-3.0 


4.50 


4.57 


+0.07 


HM, int. b 


4.44 


+0.0 


7.50 


7.52 


+0.02 


HM, flux b 


4.44 


+0.0 


7.50 


7.48 


-0.02 



The impact of 3D model atmospheres on the derived stellar metal- 
licities have been investigated for in total six Fell lines (614.9, 
623.8, 624.7, 641.7, 643.2 and 645.6 nm). In all cases the dif- 
ferent Fe II lines give the same granulation corrections to within 
0.03 dex. 

Using the Holweger-Miiller (1974) ID semi-empirical model at- 
mosphere for the solar case instead of a MARCS solar model, ei- 
ther using intensity or flux line profiles 



ture, as there still exist inherent assumptions and approxima- 
tions in the present 3D analysis, most notably the use of in- 
stantaneous mo lecu lar equilibrium and LTE radiative transfer 
(Sects. 5.5 and 5.6), besides of course the fact that no direct 



5 The results shown in Fig. [l| differ slightly from those presented 
in Asplund (2001) due to the inclusion of additional lines and simula- 
tions and use of more temporally extended 3D model atmospheres. 



comparison with observations on a star-by-star basis has been 
made. Furthermore, the issues of possible missing UV opaci- 
ties (Balachandran & Bell 1998; Bell et al. 2001) and stellar Fe 
abundances (Thevenin & Idiart 1999; King 2000, cf. Sect. |J) 
must be addressed before safe conclusions can be drawn from 
OH lines regarding the [O/Fe] behaviour in metal-poor stars. 
Therefore, we here refrain from claiming accordance between 
the UV OH and [O I] results and instead settle to point out a 
possible serious systematic error affecting recent ID analyses 
of OH lines. As a result, the case for a continuous linear in- 
crease in [O/Fe] towards lower metallicities is seemingly much 
weakened. 

Naively, one could expect that all OH molecular lines 
should be similarly affected by the low temperatures and high 
molecular number densities in the 3D model atmospheres 
of metal-poor stars. In view of this, the low [O/Fe] values 
found from OH vibrational-rotational lines in the IR from cool 
(T^ Si 5000 K) metal-poor giants and dwarfs (Balachandran 
et al. 2001; Melendez et al. 2001) appear surprising. No 3D 
hydrodynamical model atmospheres are yet available for these 
lower T c ff and/or surface gravities while the IR lines are un- 
propitiously weak for the stellar parameters of the current 
(T c ff Jb 5800 K) suites of models. Hence we are unable to con- 
firm or disprove this apparent discordance. However, we have 
verified that OH lines in the IR with similar line strengths as 
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Fig. 3. The impact of 3D model atmospheres on the [O/Fe] 
trend with metallicity for the solar (Upper panel) and turn-off 
(Lower panel) sequences. The solid line represents the ID LTE 
result of Boesgaard et al. (1999) based on only the UV OH 
lines and the King (1993) T cff -scale: [O/Fe] = -0.40-[Fe/H]. 
The dashed and dotted lines denote the typical O abundance 
corrections in a 3D LTE analysis (average for the OH 313.9 
and 316.7nm lines) compared with the ID case depending 
on the choice for the solar calibration. For the former the 
Holweger-Miiller (1974) semi-empirical solar atmosphere with 
intensity profiles has been used while for the latter a theoreti- 
cal MARCS model atmosphere together with flux profiles have 
been adopted; for all other ID models theoretical MARCS mod- 
els have been utilised. Hence, in the upper panel the granula- 
tion correction for the [Fe/H]=0.0 model is forced to disappear 
since it corresponds to the Sun 



the UV lines suffer from as severe granulation abundance cor- 
rections for the models listed in Table [l] using fake OH lines 
with strongly enhanced transition probabilities. Thus, the spec- 
tral locations of the transitions are not a solution to this conun- 
drum. Whether the differences in stellar parameters may be a 
possible resolution is discussed further in Sect. ^. 



5. Possible systematic errors affecting the 3D 
results 



Fig. 4. The impact of the 3D LTE analysis of the OH lines on 
the [O/Fe] trend in metal-poor stars when averaging the gran- 
ulation corrections for the solar and turn-off sequences for the 
two OH lines 313.9 and 316.7 nm. As in Fig. |l both the cases 
when using the Holweger-Miiller (1974) model atmosphere 
(dashed line) and the MARCS model atmosphere (dotted line) 
for the solar calibration are shown. It should be emphasized 
that the assumption of LTE in the molecule formation and ra- 
diative transfer may skew the results according to the discus- 
sion in Sect. || 



5. 1. Temperature structure in 3D model 
atmospheres 

The large granulation effects on the derived O abundances for 
OH lines in metal-poor stars hinge crucially on the existence 
of the low atmospheric temperatures encountered in the OH 
line-forming region with 3D model atmospheres. It is therefore 
natural to inquire whether the resulting temperature structures 
are indeed accurate representations of the real stellar atmo- 
spheres. Even if the here employed 3D hydrodynamical model 
atmospheres no longer rely on the same simplifying assump- 
tions as in classical ID model atmospheres, uncertainties in the 
temperature structures may still remain. A consistency check 
on the 3D atmospheres is available from similar 2D radiative- 
hydrodynamical simulations (Ludwig et al. 1999; Freytag et al. 
1999), which also produce the distinct sub-radiative equilib- 
rium temperatures. This consonance is reassuring but not too 
surprising given the resemblance in the underlying assumptions 
of the two types of simulations. 

As explained in Sect. the low atmospheric temperatures 
are natural consequences of no longer enforcing radiative equi- 
librium and instead solving explicitly the time-dependent en- 
ergy equation. This phenomenon is therefore a real physical 
effect which must be present in the photospheres of metal- 
poor stars. Nevertheless, the magnitude of the effect may have 
been overestimated in the present simulations. Asplund et al. 
(1999a) indeed cautioned that the neglect of Doppler shifts in 
the treatment of the strong spectral lines in the construction of 
the 3D model atmospheres may lead to the radiative heating 
being underestimated. Work is currently being undertaken to 
construct such further improved 3D model atmospheres. 
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An alternative avenue to proceed is to design observational 
tests to confront the predictions from 3D model atmospheres 
with. An often-used method is a detailed comparison of spec- 
tral line asymmetries and shifts (e.g. Asplund et al. 2000a,b; 
Allende Prieto et al. 2001). Unfortunately, a similar study of 
the asymmetries of the OH lines is unlikely to give crucial clues 
due to the inevitable blends in the UV. Furthermore, only lim- 
ited guidance is likely to be obtained from studies of line asym- 
metries of atomic species since they do not probe the same high 
atmospheric layers as the molecular transitions. 

A possible venture to explore is the individual O abun- 
dances obtained from a large sample of OH lines of differ- 
ent excitation potential. Since the exact 3D abundance correc- 
tions depend on the transition properties with low-excitation 
and strong lines being more influenced by the 3D model atmo- 
spheres (Table ||), different results are expected compared with 
ID model atmospheres. An examination of Table 3 in Israelian 
et al. (1998) indeed reveals that low excitation OH lines (e.g. 
312.39, 312.81, 312.83 nm, Xcxc ^ 0.2 eV) appear to systemat- 
ically give ~ 0. 1 dex higher abundance than the high excitation 
lines (e.g. 316.71, 320.39, 325.55 nm, Xc xc ^ 0.8 - 1.3eV) in 
a ID analysis, hinting to a possible problem in the ID tempera- 
ture structures at low metallicities relative to the Sun. However 
there is no clear trend with metallicity in the differences in 
the derived O abundances but the scatter is unfortunately very 
large. A proper investigation will therefore require many more 
lines for a larger stellar sample. Furthermore, a prerequisite is 
an accurate determination of the stellar parameters, which, in 
view of the current disagreement on the T c ff-scale of metal- 
poor stars, is unlikely to be settled for good in the near future. 



5.2. Effective temperature 

In order to obtain accurate absolute abundances, not only a re- 
alistic model of the stellar photosphere and a proper under- 
standing of the line formation process are necessary but also 
appropriate fundamental stellar parameters. For OH lines, the 
effective temperature is of special importance. In Sect, [j] the 
comparison between the 3D and ID predictions were carried 
out assuming that the relevant T c g should be the same for the 
two types of models. 

One may suspect that the presence of temperature inho- 
mogeneities in the continuum-forming layers should make the 
emergent flux distribution different in 3D model atmospheres 
compared with homogeneous ID models, in particular for 
metal-poor stars which often are characterized by "naked gran- 
ulation" (Nordlund & Dravins 1990; Asplund et al., in prepara- 
tion): the region of maximum horizontal temperature contrast 
reaches the visible surface whereas for example for the Sun 
these layes are hidden slightly below the photosphere (Stein 
& Nordlund 1998). This difference will be particularly man- 
ifested in UV colours (amounting to < 10% in continuum 
fluxes at 314nm, cf. Fig. [j]) while the effect is minimized at 
IR wavelengths. Since IRFM is our preferred choice for T e g- 
calibrations, we will here only investigate the impact of 3D 
models on this method. 



The IRFM is designed to compare the observed ratio of the 
total bolometric flux of a star and a monochromatic continuum 
flux at IR wavelengths with the corresponding theoretical ratio 
from model atmospheres (Blackwell & Shallis 1977), In prac- 
tice, JHK photometry normally replaces the monochromatic 
continuum flux for observational convenience. Here we follow 
the original idea by computing the spatially and temporally av- 
eraged continuum flux at 2.2 fim for the various 3D simula- 
tions and compare these predictions with corresponding ones 
for ID model atmospheres with T e g 3D and T e g^ ± 100 K 
to quantify typical corrections to T c ff-estimates based on the 
IRFM and classical ID model atmospheres. It should be noted 
that our procedure thus neglects the model atmosphere depen- 
dence in estimating the stellar fluxes outside the observed pho- 
tometric bands to obtain bolometric fluxes (Alonso et al. 1995). 
However, only a small fraction (^ 10%) of the total flux is car- 
ried at those wavelengths for the F-K stars of interest here and 
the slight inconsistency of relying on ID models for this pur- 
pose should have a marginal effect on the final T c g-calibration. 

The differences in T c s -determinations from 3D and ID 
model atmospheres are shown in Fig. 0. Clearly due to the 
small model sensitivity of the method, the use of ID models 
in available IRFM determinations does not significantly en- 
cumber the results. For the solar sequence the T c s differences 
amount to typically 20 K while the modifications for the turn- 
off sequence are completely negligible. The smaller effect for 
the hotter models is in accordance with the relatively smaller 
fraction of the total flux emitted at IR wavelengths and thus 
smaller sensitivity of the temperature structure of the adopted 
model atmospheres. In view of the typical observational un- 
certainties of « 100 K currently attached to IRFM, the errors 
in IRFM T c ft estimates are still very much dominated by the 
accuracy of the observations rather than the adopted model at- 
mospheres. 

The small influence on T e g translates to only a minor im- 
pact on the estimated O abundances. An increase in T c g by 
+100 K typically implies an increase of the derived O by 
0.2 dex for the UV OH lines (Nissen et al. 1994). Thus, a new 
Teg-calibration based on 3D model atmospheres is only ex- 
pected to introduce a 0.05 dex alteration of the inferred O 
abundances. Furthermore, there is no significant metallicity de- 
pendence in the Jeff-corrections, which could bias any deduced 
[O/Fe] trends. 

We conclude that [O/Fe] determinations will not be sub- 
stantially modified due to changes in the stellar parameters en- 
tailed by the adoption of 3D model atmospheres. 



5.3. Surface gravity 



In addition to a specification of T c ff (Sect. 5.2), knowledge of 



the stellar surface gravity, log g, is required for the spectral syn- 
thesis. In Sect. ^| the 3D and ID calculations were performed 
with identical log g, which may not be appropriate. 

The best method to determine stellar surface gravities is 
to make use of the accurate parallaxes now available from the 
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Fig. 5. The differences in T c g derived from IRFM for the solar 
(solid line) and turn-off (dashed line) sequences of models in 
the sense T^jd- T e g t ±D. The dotted line represent equality 
between the T c g determinations 

Hipparcos mission. By manipulating the familiar relationships 
g oc M/R 2 and L oc R 2 T^ S , the trigonometric gravities are: 

\ogg = 41ogT c ff+logAf/M Q +21og7r+0.4(y+SC)-10.51(l) 

where M is the stellar mass, V is the apparent visual magni- 
tude of the star corrected for extinction, BC is the bolometric 
correction and ir is the parallax (e.g. Nissen et al. 1997). The 
additional uncertainty introduced by the use of 3D model at- 
mospheres instead of classical ID models enters through the 
small systematic differences in T e g and BC, while the indirect 
effect on the derived mass through isochrone-fitting is negli- 
gible. As seen in Sect. 5.2, the difference in T^ff-calibration 
between the two types of models amount to ±20 K in the 
relevant temperature range, which corresponds to an error of 
merely ±0.01 dex for log g. No bolometric corrections have 
yet been computed with 3D model atmospheres but consid- 
ering that AT off = 500 K corresponds to ABC ~ 0.05 for 
T e g ~ 6000 K with ID model atmospheres (e.g. Alonso et al. 
1995), a very conservative limit to the differential effect of 3D 
models is ABC < 0.1. Thus, log g determinations using 3D 
model atmospheres are only expected to be < 0.05 dex dif- 
ferent from a corresponding ID analysis when relying on the 
parallax-method. We emphasize that this is not the full uncer- 
tainty with this method, which is dominated by the errors in 
the parallax measurements and the T e g -calibration, but only 
the additional uncertainty relative to the ID case when relying 
on 3D model atmospheres. In terms of derived O abundances, 
Alogg < ±0.05 corresponds to A(0/H) < =F0.02dex for the 
UV OH lines (Nissen et al. 1994). We can therefore safely con- 
clude that our assumption of identical adopted surface gravities 
in the 3D- ID comparison in Sect. |] will not impede the con- 
clusions presented therein. 

Alternative methods to derive stellar surface gravities are 
available from spectroscopy, in particular by enforcing ioniza- 
tion equilibrium or using the pressure-damped wings of strong 
lines. Unfortunately they suffer from several drawbacks which 
make them less attractive in analyses with 3D model atmo- 
spheres. In LTE the 3D abundances derived from Fe I lines 



are much lower than in classical ID analyses (Asplund et al. 
1999a), which, if correct, would lead to large modifications 
of the derived log g from ionization balance. However, the 
Fe I lines are almost certainly seriously affected by departures 
from LTE in 3D models, and thus a 3D non-LTE study of Fe 
line formation would be required, a very challenging task for 
the future. Indeed, even with ID model atmospheres depar- 
tures from LTE make the ionization gravities discrepant from 
trigonometric gravities (Allende Prieto et al. 1999b, cf. Sect. 
5.4). Similarly, the strong lines normally utilised for gravity- 
determinations are from species which can be expected to be 
affected by departures from LTE in 3D model atmospheres 
(Mg I, Ca I, Fe I), again necessitating 3D non-LTE investiga- 
tions. Furthermore, at very low metallicities ([Fe/H] —2) 
also the strongest lines become too weak to accurately probe 
the photospheric pressure structure (Fuhrmann 1998). 

5.4. Metallicity 

Stellar metallicities enter into abundance analyses both indi- 
rectly through their influence on the photospheric structure 
and directly through the use of Fe as a reference element for 
abundance ratios. While the dependence on the former is rel- 
atively weak (an error in [Fe/H] as large as 0.4 dex only im- 
plies an uncertainty in the derived O abundance from OH lines 
of 0.05 dex, Nissen et al. 1994), the latter is as important as 
deriving accurate O abundances when attempting to trace the 
evolution of [O/Fe]. This obvious fact is, however, often over- 
looked with values simply taken from the literature or esti- 
mated only from Fe I lines with no consideration for depar- 
tures from LTE. Since the recent analyses of UV OH lines by 
Israelian et al. (1998) and Boesgaard et al. (1999) have utilised 
Fe I lines, their derived [O/Fe] trends may be systematically 
overestimated (King 2000; Israelian et al. 2001), which in turn 
could influence the 3D- ID comparison presented in Sect. [|. 
Recently two investigations of non-LTE effects in Fe I line for- 
mation in metal-poor stars have been published (Gratton et al. 
1999; Thevenin & Idiart 1999) although with discomfortingly 
discordant results, which deserves further scrutiny. 

We tend to view the calculations of Gratton et al. (1999) 
with some balanced scepticism. Their incomplete treatment of 
the high-excitation levels, neglect of available quantum me- 
chanical calculations for the photo-ionization cross-sections 
(from e.g. the IRON Project, Bautista 1997) and their ex- 
ceedingly large cross-sections for inelastic collisions with H, 
all combine to ensure a result close to the LTE prediction. 
Thevenin & Idiart (1999) on the other hand adopt more realistic 
atomic input data but still suffer from the incomplete handling 
of the line-blanketing. Since the main non-LTE effect, over- 
ionization, feeds on the UV radiation field it is paramount to 
address the UV line-blocking in the calculations of the photo- 
ionization rates to avoid predicting too large departures from 
LTE for Fe I. Furthermore, improved quantum mechanical cal- 
culations for the H-collisions are urgently needed to replace the 
questionable classical recipe of Drawin (1968). 

Until improved non-LTE calculations are available, we urge 
that the determinations of stellar Fe abundances to be based 
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on Fe II lines. Fe II lines are essentially immune to departures 
from LTE (e.g. Shchukina & Trujillo Bueno 2001) and as 
clear from Table || are not particularly affected by the temper- 
ature inhomogeneities and different temperature structures in 
3D model atmospheres. At this stage, however, we can not ex- 
clude that the estimated [O/Fe] trend with 3D models (Fig. |j) 
may even need further downward adjustment due to departures 
from LTE for Fe I lines, on which the existing ID [Fe/H] esti- 
mates are based (Israelian et al. 1998, 2001; Boesgaard et al. 
1999). Naturally, investigations of departures from LTE for Fe 
in metal-poor stars should also be based on the new generation 
of 3D model atmospheres. 

5.5. Molecular equilibrium for OH formation 

For the LTE line formation calculations presented in Sect. the 
assumption of instantaneous molecular equilibrium has been 
made in the computation of the total number density of OH 
molecules at different times and locations in the 3D model 
atmospheres. According to Fig. [j], the OH line formation is 
strongly biased towards the upfiow regions, where the gas is 
very rapidly cooled from about 10000K to about 4000 K in 
a relatively thin zone around continuum optical depth unity. 
This transition occurs on a time-scale of merely a few minutes 
for the upflowing material in the Sun, which could imply that 
molecular equilibrium is not established. Additionally, pho- 
todissociation due to the non-local radiation field from deeper 
layers may cause further departures from LTE. As a result, one 
would expect that LTE may overestimate the OH content and 
therefore that the O LTE abundances may be underestimated. 

In principle it is a straight-forward exercise to compute the 
OH molecule formation and the resulting OH number densities 
as a function of time in our 3D model atmospheres by solv- 
ing a set of coupled differential equations corresponding to a 
network of chemical pathways. However, a major obstacle is 
the apparent lack of rate coefficients for the relevant reactions 
and temperatures, both experimental and theoretical. We have 
scoured various publically available databases such as UMIST 
(Le Teuff et al. 2000) in search of the necessary rate coeffi- 
cients but with little success in locating data for T J; 3000 K, 
which prevents us from performing the non-LTE chemistry cal- 
culations. 

Some guidance to the non-LTE behaviour may still be ob- 
tained from observations in the absence of detailed compu- 
tations. Uitenbroek (2000a,b) has recently concluded from a 
comparison of the observed solar CO line intensities and their 
temporal variations with calculations based on both ID hy- 
drodynamical chromospheric simulations (Carlsson & Stein 
1992, 1995, 1997) and 3D model atmospheres similar to those 
utilised here (Stein & Nordlund 1998) that the inherent as- 
sumption of instantaneous molecular equilibrium for CO may 
not be valid in the Sun. This may suggest that a similar phe- 
nomenon could also occur for OH. If so, it would probably 
be more pronounced in metal-poor stars in view of their more 
rapid and dramatic cooling in the photosphere. 

We conclude that it can not be excluded that departures 
from molecular equilibrium may influence the derived O abun- 



dances when using OH lines and that its magnitude could be 
metallicity dependent. Naturally this should be examined fur- 
ther once the necessary data becomes available. 

5.6. Local thermodynamic equilibrium for OH line 
formation 

Besides the assumption of chemical equilibrium for the OH 
molecule formation, it should be borne in mind that LTE has 
also been assumed in the OH line transfer calculations pre- 
sented in Sect. |j. Due to the vast number of relevant levels and 
transitions, detailed non-LTE calculations for the OH radiative 
transfer is formidable. In fact, non-LTE effects for molecules in 
stellar atmospheres is largely unexplored even with ID model 
atmospheres. 

We are not aware of any non-LTE studies for OH but CO 
has attracted slightly more attention. Uitenbroek (2000b) has 
recently performed detailed non-LTE calculations for the CO 
vibrational-rotational lines in the Sun, which confirm the in- 
sightful prediction by Hinkle & Lambert (1975) that the lines 
are collisionally controlled and thus that LTE is a good ap- 
proximation. Hinkle & Lambert caution on the other hand that 
electronic molecular transitions, like the UV OH lines, may 
be radiatively determined. If one approximates the UV OH 
line formation with the two-level approach with complete re- 
distribution (cf. Mihalas 1978), the line source function S\ 
will depend on the mean intensity averaged over the absorp- 
tion profile J\ = J (f)\J\d\ and Planck function B\(T) as 
S\ = (1 — e) J\ + eB\(T). Here e is a measure of the photon 
destruction probability (0 < e < 1). In the UV, J\ tend to be 
larger than B\{T) for weak lines and therefore S\ > B\{T). 
As a consequence, one would expect the OH lines to be weaker 
with scattering than in LTE, or, equivalently, that the derived O 
abundance will be underestimated in LTE. 

A detailed non-LTE calculation for OH including all the 
vibrational and rotational levels is unfortunately beyond the 
scope of the present investigation. But we have nevertheless 
attempted to estimate the non-LTE corrections to the derived 
O abundances from the UV OH lines using a two-level OH 
molecule. Although no doubt unrealistically simplistic to allow 
accurate quantitative estimates of the non-LTE effects, the ap- 
proach is still expected to yield qualitatively correct results of 
the non-LTE behaviour for different stellar parameters. For the 
purpose, version 2.2 of the statistical equilibrium code MULTI 
(Carlsson 1986) has been used after some minor modifications 
to allow treatment of OH molecular lines. The two levels cor- 
respond to the OH 313. 9nm transition with the same adopted 
parameters as for the LTE calculations presented in Sect. || 
Additionally, cross-sections for collisions with electrons and 
hydrogen must be specified. For the former the classical recipe 
of van Regemorter (1962) was adopted in the absence of more 
appropriate treatments. For H-collisions the calculations were 
performed with the formula by Drawin (1968) multiplied by a 
variable factor x. It should be noted that the Drawin formula 
was developed for atoms and it is not clear whether it is at all 
applicable for molecules. Although often used in non-LTE cal- 
culations for late-type stars (cf. discussion in Kiselman 2001), 
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Table 4. The ID non-LTE O abundance corrections when us- 
ing a two-level model-molecule to study the UV OH line for- 
mation. The last two ID models correspond to HD 140283 and 
G64-12, respectively. 
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+1.2 


0.0 
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+0.5 
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+0.25 










0.01 


+0.25 










1.0 


+0.24 








+1.0 


0.0 


+0.27 










0.01 


+0.27 










1.0 


+0.26 


6450 


4.04 


-3.0 


+1.0 


0.0 


+0.15 










0.01 


+0.15 










1.0 


+0.15 



The solar O abundance is here assumed to be logeo = 8.90 
The factor which is multiplied to the Drawin (1968) recipe for the 
collisional cross-sections with H 

The non-LTE O abundance correction defined as Alogeo 

LT 



logeg LTE 



log^ TE 



evidence is now mounting that the simple-minded approach by 
Drawin severely overestimates the importance of H-collisions 
by at least three orders of magnitude for atomic transitions 
(Fleck et al. 1991; Belyayev et al. 1999). The inclusions of vi- 
brational and rotational sub-levels and line-blanketing in the 
UV radiation field are also expected to diminish the non-LTE 
effects due to stronger collisional quenching and decreased 
J\/B\(T) ratios in the line-forming region. 

Non-LTE calculations have been performed for four ID 
model atmospheres with the solar T c s and log g but different 
metallicities ([Fe/H] = 0.0, -1.0, -2.0, and -3.0), as well 
as for selected metal-poor stars (e.g. HD 140283 and G64-12) 
for different O abundances to estimate the influence of the line 
strength. In all cases, the effect of H collisions were investi- 
gated by computing the abundance corrections for three differ- 
ent x: 0.0, 0.01 and 1.0. The ID models are not identical to 
those adopted for the solar sequence in the 3D- ID LTE com- 
parison in order to strictly isolate the metallicity dependence of 
the non-LTE effects. 




313.875 313.880 313.885 313.890 313.895 
Wavelength [nm] 

Fig. 6. The ID OH313.9nm line in a metal-poor ([Fe/H] = 
—3.0) Sun in LTE (solid lines), and in non-LTE with x = 1.0 
(dashed lines) and x = 0.0 (dotted lines), see text for de- 
tails. Clearly the inclusion of collisions by H according to the 
Drawin (1968) formula (x ~ 1.0) has a negligible impact on 
the resulting line strength 



A summary of the results in terms of abundance corrections 
is presented in Table [|. As expected, in all cases the assump- 
tion of LTE makes the lines stronger, and therefore that LTE 
analyses underestimate the O abundances. The difference in 
terms of line strengths is quite large (Fig. ||) and therefore rel- 
atively large non-LTE corrections are obtained, ~ 0.2 dex. The 
predicted ID non-LTE effects are almost independent on the 
stellar parameters. Fortunately, the magnitude of the non-LTE 
corrections is only marginally dependent on the treatment of 
the H collisions provided the Drawin (1968) formula does not 
greatly underestimate the collisional cross-sections {x > 10). 
Even with x = 1.0 the predicted non-LTE corrections differ 
only by < 0.03 dex compared to the case when neglecting the 
H collisions completely. In all cases, the same results are ob- 
tained with x = 0.01 and x — 0.0. We emphasize though that 
it is still unclear whether the Drawin recipe can be applied also 
to molecules. But at this stage there are no indications that the 
treatment of H collisions play any significant role in the pre- 
dicted non-LTE corrections. 

The departures from LTE in this two-level approach 
are purely due to scattering in the line, in agreement with 
the prediction by Hinkle & Lambert (1975) (Fig. 0). Since 
J\/B\(T) > 1 at the relevant atmospheric layers for the stud- 
ied model atmospheres, the line source function S\ exceeds the 
LTE value. This effect is more pronounced in the metal-poor 
models but its influence on the abundance corrections is partly 
balanced by the weaker lines at lower metallicities. Since the 
vast majority of the OH molecules will be in the ground elec- 
tronic state and the vibrational and rotational levels within elec- 
tronic states are closely coupled through collisions, the depar- 
ture coefficient j3 = Ni/Nf' TE for the lower level will be very 
close to one and thus no non-LTE effects due to differences in 
line opacity emerge, as evident from Fig. ^| 

In spite of the significant ID non-LTE abundance correc- 
tions, we find no evidence for a pronounced steepening of the 
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Fig. 7. The ratio of the line source function with the Planck 
function S\/B\(T) for the two cases x = 1.0 (dashed lines) 
and x — 0.0 (dotted lines) in a metal-poor ([Fe/H] = —3.0) 
Sun, see text for details. Also shown are the departure coeffi- 
cients (3 = iVi/iV 1 LTE for the lower level of the OH transition 
for x — 1.0 and x — 0.0 (solid lines), emphasizing that the 
non-LTE effect is purely due to scattering and not changes in 
the line opacity 



[O/Fe] trend. However, we caution that this conclusion can not 
automatically be extrapolated to the 3D case, in which the de- 
partures from LTE may be more severe in metal-poor stars 
since the steep temperature gradients may be more prone to 
scattering effects than in ID model atmospheres. Needless to 
say, an investigation of the non-LTE behaviour in 3D model 
atmospheres as a function of metallicity has very high priority. 

5.7. Summary of uncertainties: 3D or ID model 
atmospheres? 

Considering that the results presented here are among the first 
investigations of the impact of the new generation of 3D hy- 
drodynamical model atmospheres on stellar spectroscopy, it is 
therefore in order to ask if indeed these models are more realis- 
tic than previously used classical ID model atmospheres. Sofar 
the predictions from the ab-initio 3D models have been very 
successfully confronted with detailed observational constraints 
for in particular the Sun. These comparisons include such dis- 
parate tests as granulation topology and flow field (Stein & 
Nordlund 1998), helioseismology (Rosenthal et al. 1999), in- 
tensity brightness contrast (Stein & Nordlund 1998; Asplund 
et al. 2000a), flux distribution and limb-darkening (Asplund 
et al. 1999b) and spectral line shapes, shifts and asymmetries 
(Asplund et al. 2000b). No doubt the current surface convection 
simulations for the Sun have a very high degree of realism. In 
sharp contrast theoretical ID model atmospheres fail in regards 
to most, if not all, of the above-mentioned tests. Recently sim- 
ilar 3D models have been used for studies of line asymmetries 
in Procyon (Allende Prieto et al. 2001) and the metal-poor halo 
star HD 140283 (Allende Prieto et al. 1999a; Asplund et al., in 
preparation) with very satisfactory outcomes. 



In view of the detailed discussion given in Sect. |5| on pos- 
sible remaining systematic errors in the 3D analysis, the reader 
may get the impression that the here presented results are rather 
uncertain. It may therefore be in order to point out that all 
of these possible effects also apply equally well to any study 
based on ID model atmospheres, besides the errors introduced 
by the assumption of hydrostatic equilibrium and by treating 
convection through the mixing length theory (or a close rel- 
ative thereof), which is known to be a poor representation of 
stellar convection. However, at this stage it is probably prema- 
ture to conclude that 3D model atmospheres are indeed supe- 
rior to classical ID models. It is therefore of utmost importance 
now to carry out the same arsenal of tests which previously has 
been undertaken with ID model atmospheres (flux distribution, 
limb-darkening, colours, H-lines etc) as well as additional ones 
now possible (line asymmetries and shifts, asteroseismology), 
in particular for metal-poor stars. This is even more true in light 
of the fact that ID model atmospheres often fail the very same 
tests. 

6. Comparing previous OH, Oi and [Oi] results 

Even if the present article does not deal directly with observa- 
tions, in this section we will nevertheless discuss some of the 
recent analyses of various O diagnostics as it has bearing on our 
findings and may give clues to the existence or not of the large 
granulation corrections for OH in metal-poor stars described in 
Sect. £§ 

A major argument for the monotonic linear trend in [O/Fe] 
with metallicity derived from the UV OH lines comes from the 
claimed good agreement with the O I triplet results (Israelian 
et al. 1998, 2001; Boesgaard et al. 1999). Since much smaller 
granulation corrections are expected for the triplet than for OH 
(Asplund 2001) this would seem to contradict the findings in 
Sect. However, the consonance between the OH and O I re- 
sults is not as unambigous when examining some of the pub- 
lished analyses in detail. Fig. [| shows the difference between 
the OH-based and 1-based abundances of Boesgaard et al. 
(1999) on the King (1993) T cff -scale; with the lower Carney 
(1983) scale the correlation is slightly less pronounced but 
clearly present (slope = —0.13 instead of —0.18). The abso- 
lute value for the abundance differences is here less important 
than the existence of the metallicity-trend due to the uncertainty 
introduced by the choice of T ff-scale and the neglect of non- 
LTE effects for O I. This divergent behaviour is very close to 
the expected according to the results of Sect. ^ and Asplund 
(2001), which we interpret as a qualitative argument for the 
aptness of the 3D calculations. In view of this, the good gen- 
eral agreement found by Israelian et al. (1998) for nine stars 
using the Oi equivalent widths of Tomkin et al. (1992) pose 
a perplexing problem. Clearly a larger stellar sample with si- 
multaneous analyses of the OH and O I lines would be very 
helpful in this context. One should also explore possible differ- 
ences in for example C abundances derived from C I and CH 
lines at low metallicities. Preliminary calculations reveal sig- 
nificant but smaller differences between ID and 3D analyses 
of CH lines compared with for OH lines (« 0.3 dex instead of 
« 0.6 dex for OH at [Fe/H] = -3.0). 
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Fig. 8. The difference between the OH and 1-based ID abun- 
dances of Boesgaard et al. ( 1 999) on the King ( 1 993) -scale, 
which reveals a pronounced metallicity-dependent trend (slope 
—0.18). This behaviour is in qualitative agreement with the 3D 
calculations presented in Sect. |] 



The conventional [O/Fe] plateau indicated by the [O I] line 
has recently received support from the IR OH vibrational- 
rotational lines (Balachandran et al. 2001; Melendez et al. 
2001). As discussed in Sect. this could be interpreted as the 
here presented granulation corrections for the UV OH lines in 
metal-poor stars are overestimated, since for the same star the 
IR and UV lines should to first order be affected similarly (with 
the possi ble e xception of non-LTE effects in the line formation, 
cf. Sect. 5.6). A more likely explanation in our opinion, how- 
ever, is to be found in the differences in stellar parameters, in 
particular T e g, for the sofar investigated UV and IR targets. 
As evident from Table || and Fig. || the magnitude of the OH 
granulation corrections depend on T c g in the sense that larger 
effects are present for higher temperatures. Although danger- 
ous to extrapolate, we speculate that significantly smaller OH 
abundance corrections than found here will be obtained for 
T e s ^ 5000 K, which is typical for the metal-poor stars with 
detected IR OH features. The physical reason for the depen- 
dence on T c s is the strength of the coupling between the gas 
and radiation field: for lower T c g, more and stronger spectral 
lines are available which contribute additional radiative heat- 
ing, keeping the gas temperature closer to the radiative equilib- 
rium value (Asplund et al., in preparation). No 3D model at- 
mospheres are yet available for these lower T e g which prevent 
verification of our hypothesis but such convection simulations 
are currently being constructed. 

Due to the inconspicous nature of the [Ol] line (Nissen 
et al. 2001), no direct comparison with the UV OH results 
has sofar been possible for metal-poor dwarfs with [Fe/H] < 
—2. Recently, Fulbright & Kraft (1999) have considered the 
[Ol] line in the two metal-poor ([Fe/H] < —2.1) subgiants 
BD+23°3130 and BD+37°1458 (T cS = 5130K and 5260 K, 
respectively), which have been re-analysed by Israelian et al. 
(2001). The UV OH lines suggest [O/Fe] = 0.60 and 0.50, re- 
spectively, in the two stars, i.e. in significantly better agreement 



with the traditional [O/Fe] plateau compared with the linear 
trend seemingly implied by the OH lines in the hotter stars. As 
for the IR analyses, we conjecture that at T c g ~ 5200 K signifi- 
cantly smaller granulation corrections than those given in Sect. 
[| will be deduced, leaving the plateau-like [O/Fe] ratios largely 
unchallenged. Excellent agreement between [O I] and OH is 
found for both stars: the forbidden line gives [O/Fe] = 0.62 
and 0.52, respectively. The latter values come from our own 
analysis using ID MARCS model atmospheres with the stel- 
lar parameters given in Israelian et al. (2001). The equivalent 
widths for [O I] and Fe II were taken from Fulbright & Kraft ex- 
cept for [O I] in BD+23°3130 for which we relied on the VLT- 
measurement (W x = 0.15pm) by Cayrel (2001). This [O/Fe] 
estimate for BD+23°3130 is for unknown reason significantly 
lower than reported in Israelian et al. (2001), [O/Fe] = 0.82, 
while in perfect agreement with the findings of Balachandran 
(2001, private communication) using Kurucz (1993) model at- 
mospheres but otherwise identical input. 



7. Concluding remarks 

The low atmospheric temperatures encountered in the new gen- 
eration of 3D hydrodynamical model atmospheres (Asplund et 
al. 1999a, 2000a,b, Allende Prieto et al. 2001; Asplund et al., 
in preparation) compared with classical ID hydrostatic model 
atmospheres for metal-poor stars have a profound impact on 
the OH line formation, as described in Sect. Q As a conse- 
quence, a possible severe systematic error in recent ID LTE 
analyses of OH lines in metal-poor stars (Israelian et al. 1998, 
2001; Boesgaard et al. 1999) has been exposed, making the 
case for a monotonic increase in [O/Fe] towards lower metal- 
licities less convincing. Taken at face value, our 3D LTE OH 
analysis results in [O/Fe] values in rough agreement with the 
conventional [O/Fe] plateau for [Fe/H] ^ — 1 indicated by the 
[Ol] lines (Fig. Q). We emphasize though that this apparant 
concordance should not be taken too literally in view of the 
preliminary nature of our 3D calculations and that no star-by- 
star comparison has been made with observations. By investi- 
gating possible systematic errors for the 3D LTE results it has 
been found, however, that the conclusion of large granulation 
effects on the OH lines appears reasonably robust. For example, 
r c ff-calibrations using 3D model atmospheres should not dif- 
fer significantly (AT ff <> 20 K) from previous ID calibrations 
provided they are based on IRFM; the same is not necessarily 
true for alternative calibrations using Balmer lines or colours. 
Similarly, estimates of log g and [Fe/H] should remain essen- 
tially unaltered when relying on Hipparcos parallaxes and Fe II 
lines. The major remaining uncertainties in the 3D analysis ap- 
pear to be the assumptions of LTE for the molecular equilib- 
rium and in the line formation. Although not yet investigated 
for 3D model atmospheres, it is possible that such departures 
from LTE may indeed steepen the [O/Fe] trend once again but 
unlikely as much as the original ID LTE case. On the other 
hand, departures from LTE for Fe I may diminish the slope fur- 
ther. A final verdict on this issue must therefore await detailed 
and improved non-LTE calculations. But we note in the mean- 
time that even in the presence of possible non-LTE effects for 
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OH and Fe I, [O/Fe] ratios will be less affected since both the 
O and Fe LTE abundances will tend to be underestimated. 

The purpose of the present paper has not been to advocate 
a specific [O/Fe] trend with metallicity since our investigation 
has been limited to a differential 3D- ID comparison without in- 
volving observational confrontation. Nevertheless, our findings 
will likely fuel the long-standing debate on the O abundances 
in metal-poor stars. As already stated several times, it would 
be premature to conclude from our analysis that the mono- 
tonic linear trend in [O/Fe] claimed by Israelian et al. (1998, 
2001) and Boesgaard et al. (1999) must now be abandoned, 
even if the arguments for continuously increasing [O/Fe] ratios 
towards lower metallicities appear much weaker. It should be 
remembered that many, albeit not all, studies of the O I triplet 
at 777 nm in dwarfs and subgiants find [O/Fe] values system- 
atically higher than those given by the [O I] line in metal-poor 
giants and subgiants. The triplet is unfortunately sensitive to 
the adopted Toff-scale and departures from LTE, which de- 
serve very careful treatment. For the moment, the most reliable 
[O/Fe] ratios still appear to come from the [O I] line, provided 
that very high S/N spectra are utilised and the metallicities are 
estimated from Fe II lines (Nissen et al. 2001; Lambert 2001). 

To ignore the existing systematic errors in traditional ID 
LTE analyses of the UV OH lines in the hope that departures 
from LTE in the molecule formation and radiative transfer will 
conspire to exactly compensate the effects of the low atmo- 
spheric temperatures in 3D is certainly fraught with danger and 
furthermore most likely misleading. As demonstrated above, in 
order to derive reliable oxygen abundances in very metal-poor 
stars, unfortunately there is no easier escape route but to per- 
form time-dependent non-LTE calculations in 3D model atmo- 
spheres for OH. Fortunately, recent improvements in the anal- 
ysis of stellar spectra ensure such studies in fact being tractable 
tasks for the near future, which should help remove many of 
the lingering uncertainties in the derived oxygen abundances 
and perhaps finally settle the long-standing debate on [O/Fe] in 
metal-poor stars. 
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